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Abstract 


It is known that smooth solutions to the non-isentropic Navier-Stokes equations 
without heat-conductivity may lose their regularities in finite time in the presence 
of vacuum. However, in spite of the recent progress on such blowup phenomenon, 
it remain to give a possible blowup mechanism. In this paper, we present a simple 
continuation principle for such system, which asserts that the concentration of the 
density or the temperature occurs in finite time for a large class of smooth initial 
data, which is responsible for the breakdown of classical solutions. It also give an 
affirmative answer to a strong version of conjecture proposed by J.Nash in 1950s. 
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1 Introduction 

In this paper, we consider the system of partial differential equations for the three- 
dimensional compressible, and non-isentropic Navier-Stokes equations in the Eulerian 
coordinates 


p t + di v(pu) = 0, 

< ( pu)t + di v(pu <S) u) — pAu — (p + A)Vdivu + VP = 0, 


( 1 . 1 ) 


^c v [(p8)t + div(pw0)] — nA9 + Pdivu = 2p\T)(u)\ 2 + A(divu) 2 , 

where t > 0 is time, x € U C M 3 is the spatial coordinate, and p,u = (ui,U 2 ,u ^) tT , 9, 
P = Rp9 (R > 0), represent respectively the fluid density, velocity, absolute tempera¬ 
ture, pressure; D(u) is the deformation tensor given by 



The constant viscosity coefficients p and A satisfy the physical restrictions 


p > 0, 2p + 3A > 0. 


( 1 . 2 ) 
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Positive constants c„, k. and v are respectively the heat capacity, the ratio of the heat 
conductivity coefficient over the heat capacity. 

The compressible Navier-Stokes system m consists of a set of equations describing 
compressible viscous heat-conducting flows. Indeed, the equations CEZED i , ( 0 ) 2 , and 
(0)3 respectively describe the conservation of mass, momentum, and energy. 

There is a considerable body of literature on the multi-dimensional compressible 
Navier-Stokes system O) by physicists and mathematicians because of its physical 
importance, complexity, rich phenomena, and mathematical challenges; see [2j|2j,[5l 
[9l[l3l[l5HIZll20l[2ll[25] and the references cited therein. However, many physically 
important and mathematically fundamental problems are still open due to the lack of 
smoothing mechanism and the strong nonlinearity. For example, although the local 
strong solutions to the compressible Navier-Stokes system m for general initial data 
with nonnegative density were respectively obtained by [2], whether the unique local 
strong solution can exist globally in time is an outstanding challenging open problem 
in contrast to the isentropic case [13]. 

In the presence of vacuum, as pointed out by Xin [25], non-isentropic Navier-Stokes 
equations without heat-conductivity will develop finite time singularity, see also refer¬ 
ence [1]. Indeed, very recently, Xin-Yan [26] further proved that any classical solutions 
of viscous compressible fluids with or without heat conduction will blow up in finite 
time, as long as the initial data has an isolated mass group. Their results hold for 
the whole space and bounded domains, yet the blowup mechanism is not clarified. It 
is the main purpose if this paper to resolve this key issue. Theorem 11.11 reveals that 
the concentration of the density or the temperature must be responsible for the loss of 
regularity in finite time. 

Although vacuum will lead to breakdown of smooth solutions in finite time , it is also 
important to study the mechanism of blowup and structure of possible singularities of 
general strong (or smooth) solutions to the compressible Navier-Stokes system. 

The pioneering work can be traced to Serrin’s criterion m on the Leray-Hopf weak 
solutions to the three-dimensional incompressible Navier-Stokes equations, which can 
be stated that if a weak solution u satisfies 

u e L s ( 0, T; L r ), - + - < 1, 3 < r < oo, (1.3) 

S T 

then it is regular. 

Recently, Huang-Li-Xin nn extended the Serrin’s criterion (11.31) to the barotropic 
compressible Navier-Stokes equations and showed that if T* < oo is the maximal time 
of existence of a strong (or classical) solution (p, u ), then 

{\\p\\l°°(0,T-,L°°) + IMIl s ( 0,T;L r )) = °0> (1-4) 

with r and s as in (11.31) . For more information on the blowup criteria of barotropic 
compressible flow, we refer to mmmmm and the references therein. 

When it comes to the full compressible Navier-Stokes system m, the problem 
is more complicated. In [18], Nash proposed a conjecture on the possible blowup of 
compressible heat-conductive flows. He wrote ’’This should give existence, smoothness, 
and continuation (in time) of flows, conditional on the non-appearance of certain gross 
type of singularity, such as infinities of temperature or density.” 

Under the condition that 


A < 7p, 


(1.5) 
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Fan-Jiang-Ou [3] obtained the following blowup criterion 

T—(11 ^ 11 1,00 + ll^lllRO,T;L°°)) = °°- 

Recently, under just the physical restrictions (11.21) . Huang-Li ]8] and Huang-Li-Xin [12] 
established the following blowup criterion: 

(II^III/ 2 (0,T;L°°) + II®(' u )IIl 1 (0,T;L°°)) = °°> 

where D(u) is the deformation tensor. 

Later, in the absence of vacuum, Sun-Wang-Zhang [24] established the following 
blowup criterion for bounded domains with positive heat-conductivity k > 0 that 

T ^*(ll(/h ~)IIl°°(0,T;L°°) + II#IIl°°(0,T;L°°)) = OO, (1.6) 

provided that (11.21) and (11.51) hold true. As a consequence, Nash’s conjecture is partially 
verified as [24] can’t rule out the possibility of appearance of vacuum. 

Recently, for k > 0, we [10] establish a blowup criterion allowing initial vacuum, 
which is independent of temperature, as follows 

T lim t (||p|| i( x>( 0 ,T;Lo°) + IM|l/-(0,T;i>)) = OO, (1.7) 

where r, s satisfy m- 

As a matter of fact, the blowup criterion (11.71) further implies 

t ^JIIpI|l°°(0,T;L°°) + I|0||l°°(O,T;L°°)) = (l- 8 ) 

as long as holds true. This makes Nash’s conjecture as an immediately corollary 
for positive heat-conductivity flows. Our main goal in this paper is to give an affirmative 
answer to a strong version of Nash’s conjecture without heat-conduction. 

We will assume that k = 0, and without loss of generality, take c v = R = 1. The 
system (ED is reduced to 

Pt + di v(pu) = 0, 

< (pu)t + div(pu <S> u) — pAu — (p + A)Vdivu + VP = 0, (1.9) 

Pt + di v(Pu) + Pdivu = 2fi\D(u)\ 2 + A(divu) 2 . 

The system (11.91) is supplemented with the following initial conditions: 

(p,u,P)(x,0) = (pq,uq,Pq)(x), x€R 3 , (1.10) 

(p, a, P) satisfies the far field condition: 

(p, u, P)(x, t) —» (0, 0,0) as |x| —>• oo; (1.11) 

To state the main result, we will use the following notations and conventions. 
Notations. For 1 < p < oo and integer k > 0, the standard homogeneous and 
inhomogeneous Sobolev spaces in R 3 are denoted by: 

(lp = L p (R 3 ), W k ’P = W k ’P(R 3 ), D k ’P = {u € Ll c (R 3 ) \ V k u € L p ) , 

| pi = {w€L 6 |Vu€L 2 }, H ( ) = l 2 nnj, H k = W k > 2 
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Denote by 


/ 


ft + u- V/, 



fdx. 


The strong solutions to the Cauchy problem (11.9I) - (1 1 .111) are defined as follows. 

Definition 1.1 (Strong Solutions). (p,u,P) is called a strong solution to (11.91) in 
R 3 x (0, T), if for some qo > 3, 


\p,p)> 0, ( P ,P)eC{[0,T\-,L 1 nH 1 nw 1 M), ( Pt ,P t ) e C([0,T}-L q °), 

< u G C([0, T]; Dl n D 2 ' 2 ) n L 2 { 0, T; D 2 ’ qo ), 
u t € T 2 (0, T; D 1,2 ), Jpu t € L°°(0,T;L 2 ), 

and (p, a, P) satisfies both (jl.lj) almost everywhere in R 3 x (0, T ) and (11.101) almost 
everywhere in R 3 . 

Then the main result in this paper can be stated as follows: 

Theorem 1.1. For constant q £ (3,6], assume that (po > 0,no,Po > 0) satisfies 

{ Po ,Po)eL 1 nH 1 nw 1 ’ q , M 0 6 fijnfi 2 ' 2 , Po \u 0 \ 2 eL\ (1.12) 


and the compatibility conditions: 

- pAu 0 - (p + A)Vdivn 0 + VP 0 = y/pog, (1-13) 

with g € L 2 {R 3 ). Let ( p,u,P) be the strong solution to the compressible Navier-Stokes 
system m in R 3 . If T* < 00 is the maximal time of existence, then 


provided 


^™*(\\p\\l°°(o,t-,l°°) + II$IIl°°(o,t ; l°°)) — °°- 
p > 4A. 


(1.14) 

(1.15) 


A few remarks are in order: 

Remark 1.1. Under the conditions of Theorem, 11.11 the local existence of the strong 
solutions was guaranteed in m ■ Thus, the assumption T* makes sense. 

Remark 1.2. The main contribution of Theorem, \1.1\ asserts that Nash’s conjecture 
even holds for zero heat-conductivity flows. In the case k = 0, the formation of sin¬ 
gularity is only due to the concentration of either the density or temperature. In this 
sense, we give an affirmative answer to a strong version of Nash’s conjecture. 

Remark 1.3. Condition is only used in obtaining estimates in Lemma \3.2[ 

Remark 1.4. It’s easy to prove a same continuation principle for two-dimensional 
problem without any restrictions on p,X. Since the proof is analogous and simpler, we 
omit it for simplicity. 

We may also investigate the following different boundary conditions. 

(1) D = R 3 and constants p, P > 0, ( p,u,P ) satisfies the far field condition either 
vacuum or non-vacuum: 


(p, u, P)(x, t) —>• (p, 0, P) as |x| —>• oo; 


(1.16) 
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and initial condition 

(po-p,Po-p) « 0 €^nD 2 - 2 , P0 \u 0 \ 2 e l\ (1.17) 

(2) is a bounded smooth domain. 

u = 0 on <90. (1.18) 

For constant g £ (3,6], assume that (po > 0,uo,Pq > 0) u satisfies 

(p 0 ,P 0 ) € ^(fijn^fl), u 0 £ FT 0 1 (O)ni7 2 (O). (1.19) 

(3) O = T 3 = R 3 /Z 3 and constant q € (3,6], assume that (po > 0,uq,Pq > 0) satisfy 
(po,Po) € L 1 (T 3 )niy 1 ’ 9 "(T 3 ), uq £ H 2 (T 3 ). (1.20) 


Our next theorem asserts that Nash’s conjecture also holds for different boundary 
conditions. 

Theorem 1.2. Let (p, u, P) be the strong solution to the full compressible Navier-Stokes 
system TO together with 

(p,u,P)(x,0) = (po,u 0 ,P 0 ), x£Q, (1.21) 

and 

1. (po,uo,Po) satisfy (II. 161) - (II. 1711 for O = R 3 , 

2. ( po,uo,Po ) satisfy (II . 18|1 - (11.191) for bounded domain it, 

3. (po,uo,Po) satisfy (ll.20p for periodic domain Ll, 
and the compatibility conditions: 

- pAu 0 - (p + A)Vdmx 0 + VP 0 = y/po9, (1-22) 

with g £ L 2 (Q). If T* < oo is the maximal time of existence, then 


II PlU°°(0,T;L°°) + II ^ II L°°(0,T ;L°°)) — OO- (1-23) 

provided 

p > 4A. (1.24) 

We now make some comments on the analysis of this paper. 

Let (p, u, 6) be a strong solution described in Theorem ll.il Suppose that (11.231) were 
false, that is, 


r l™ t (llp|U oc (0,T;L°°) + I|0 ||l°°(O,T;L°°)) — M 0 < +°°- (1-25) 

One needs to show that 


SU P (II(p 5 P)lli? 1 nVF 1 -8 + l|V«||.ffi) < C. 

0 <t<T* 


(1.26) 
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Higher order derivatives estimates for above quantities then follow easily from above 
regularities. 

Let’s say a few words on the regularity criterion (11.231) . In the absence of heat- 
conductivity, the equation for the temperature changes its form from parabolic to hy¬ 
perbolic type, thus resulting in loss of regularity benefiting from the smooth effect of 
heat dissipation. But fortunately, it enjoys a same nonlinear structure as the density 
equation. Since the methods in all previous works [?,ES] depend crucially on Hoff’s a 
priori estimates. The main point is how to avoid terms involving density gradient in 
calculations. It turns out to be possible to treat the terms arising from pressure gradi¬ 
ent. Some new ideas are needed to recover all the a priori estimates, that is, instead of 
the temperature 6 and pressure P, we treat the total energy E = \\u\ 2 + 9. 

Finally, the a priori estimates on both the L^L^-norm of the density gradient and 
pressure gradient along with the L\ L£°-norm of the velocity gradient can be obtained si¬ 
multaneously by solving a logarithm Gronwall inequality based on a logarithm estimate 
(see Lemma 12.41) and the a priori estimates we have just derived. 

The rest of the paper is organized as follows: In the next section, we collect some 
elementary facts and inequalities that will be needed later. The main result, Theorem 
11.11 is proved in Section [3j 


2 Preliminaries 

In this section, we recall some known facts and elementary inequalities that will be 
used later. 

First, the following existence and uniqueness of local strong solutions when the initial 
density may not be positive and may vanish in an open set can be found in [2]. 

Lemma 2.1. Assume that the initial data (po > 0,ito,Po > 0) satisfy (|1.12D -(??). 
Then there exists a positive time Tj G (0, oo) and a unique strong solution (p, u, P ) to 
the Cauchy problem (11.11) (11.101) on R 3 x (0, Ti]. 

Next, the following well-known Sobolev inequality will be used later frequently (see 

m)- 

Lemma 2.2. For p G (l,oo) and q G (3, oo), there exists a generic constant C > 0, 
which depends only on p, q such that for f G D\ and g G L p D D 1,q , we have 

II/IIlo < C\\Vf\\ L 2 , \\g\\ L oo < CWqWlp + C\\Vg\\ Lq . (2.1) 

The following logarithm estimate from mi will be used to estimate ||Vu||l°° and 

ll(Vp, VP)|| L 2 nL9 . 

Lemma 2.3. For 3 < q < oo, there is a constant C(q ) such that the following estimate 
holds for all Vu G L 2 D D 1,q , 

||Vu||l°° < C (||divtt||x,oo + ||V x u\\ L ™) log(e + \\V 2 u\\ L i) 

+ C , ||Vu|| L 2 +C'. 

Finally, we consider the following Lame system 


pAv(x) — (p + A)Vdivu(x) = f(x), x G H, 


(2.3) 
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where v = (ui, u 2 , u 3 ), / = (/i, / 2 , /s), and g,X satisfy dH2D - 
Assume that 0 is a bounded smooth domain and 

v = 0 on <90. (2.4) 

The following logarithm estimate for the Lame system (12.31) . which can be found 
in [10] will be used to estimate ||Vu|| L oo and \\Vp\\ L 2 nLq . 

Lemma 2.4. Let g, A satisfy m- Assume that f = divg where g = (gkj) 3 x 3 with 
gkj € L 2 n W 1,q for k, j = 1, • • • , 3, r € (1, oo), and q e (3, oo). Then the Lame system 
m together with m has a unique solution v € Hq nW 2,q , and there exists a generic 
positive constant C depending only on g, A, q, r and 0 such that 

||Vu|| L r < C\\g\\ L r, (2.5) 

and 

IIVuIIloo < C (1 + ln(e + ||)||£/||l°° + HpIIz/-) ■ (2.6) 


3 Proof of Theorem 11.1 . 

Throughout the rest of the section, C will denote a generic constant depending only 
on p 0 ,u 0 ,6o,T*,M 0 ,\,g. 

We start with the standard energy estimate 

Lemma 3.1. 

dx = J QpoKI 2 + Pq) dx = E 0 . (3.1) 

Next, a high energy estimate holds under the condition (11.251) . 

Lemma 3.2. Under the condition ([1.251) . as long as g > 4A, it holds that 

p(\u\ 2 + \u\ 6 )dx + J J |Vit| 2 (l + M 2 + \u\ A )dxdt < C (3.2) 

Proof. It follows from (11.251) that 


/G p| " |2+p ) 


IIpIIlidl- < C. (3.3) 

Multiplying ( 1 . 1)2 by q\u\ q ~ 2 u , and integrating over 0, one obtains by using lemma 2.1 
that 


d_ 

dt 


j p\u\ q dx + J(q\u\ q 2 [g\Vu\ 2 + (A + /x)(divu) 2 + g(q — 2)|V|u|| 2 ] 
+ + p)(V|u| 9-2 ) • ttdivu)dx 

= q J Pdiv(|u| q_2 u)dx 

< C J ps \u\ q ~ 2 \Vu\dx 

/ l<rW* + c<.)/*!•-»& ^ 

J \u\ q ~ 2 \Vu\ 2 dx + C(e)(J p\u\ q dx)~. 


< e 

< £ 


(3.4) 











Noting that |V|it|| < |Vrt|, one gets that 

q\u\ q ~ 2 [p\S7u\ 2 + (A + /r)(divu) 2 + /i(q - 2)|V|u|| 2 ] + q(A + //)(V|u| q - 2 ) • udivu 

> q\u\ q ^ 2 [p\\/u\ 2 + (A + /i)(divu) 2 + p (q — 2)|V|u|| 2 
— (A + n)(q — 2)|V|it|| • |divu|] 

= q\u\ q ~ 2 [p\S7 u\ 2 + (A + /i)(|divu| - |v|u||) 2 ] (3.5) 

+ q\u\ q ~ 2 \ji(q - 2) - ^(A + n){q - 2) 2 ]|v|u|| 2 

> q\u\ q ~ 2 [^(q - 1) - ^(A + p){q - 2) 2 ]|V|«|| 2 

Consequently, recall q = 6 and p > 4A, the left-hand side of (13.51) is greater than 

6(/r-4A)M 4 |VM| 2 . (3.6) 


Hence, taking e small enough in (13.41) and Gronwall’s inequality implies 

ft 


J p\u\ 6 dx + J J \\7u\ 2 \u\ A dxdt < C. 
Taking <7 = 4 again, one can also prove 

J p\u\ A dx + J J \Vu\ 2 \u\ 2 dxdt < C. 

This completes the proof. 

Also, 6 is always non-negative before blowup time T*. 

Lemma 3.3. As long as do > 0, it holds that 

inf 0(x, t) > 0. 
(x,t)Git 3 x(0,T*) 

Proof. The equation for P can be rewritten as 

Pt + div(Pu) + Pdivu = F > 0. 


(3.7) 

(3.8) 

□ 


(3.9) 


(3.10) 


Since the solution is smooth , we can always define particle path before blowup time. 


^ X (x,t) = u(X(x,t),t), 


X(x, 0) = x. 

Consequently, along particle path, one has 

d 


(3.11) 


dt 


P(X (x,t),t) = —2Pdivu + F 


implies 


P(X(x, t), t) = exp(—2 / divuds ) 


P( 0) + f exp(2 f divudr)Fds 
Jo Jo 


Hence, 8 > 0 follows immediately from (13.131) . 


(3.12) 

> 0. (3.13) 

□ 
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Before proving Theorem 11.11 we state some a priori estimates under the condition 

CL25D - 

Let E be the specific energy defined by 


E = 6 + 


\u\ 


(3.14) 


Let G, co be the effective viscous flux, vorticity respectively given by 

G = (2 p + A)divu — P, co = curlu. 

Then, the momentum equations can be rewritten as 

pii = VG — /iV x oo. 

Then, we derive the following crucial estimate on the L°°(0, T; L 2 )-norm of Vu 
Lemma 3.4. Under the condition (|1.25l) . it holds that for 0 < T < T*, 

I \S/u\ 2 dx + [ 

0 <t<T 


(3.15) 

(3.16) 


sup f \'Vu\ 2 dx + f f p\ii\ 2 dxdt<C. 
) <t<T J Jo J 


(3.17) 


Proof. First, multiplying (II.IIL by u t and integrating the resulting equation over 0 
show that 

J (/r|Vu| 2 + (p + A)(divu) 2 ) dx + J p\u\ 2 dx 
= J pu- (u- V)udx + J Pdivutdx 

< - j p\u\ 2 dx + C J p\u\ 2 \Vu\ 2 dx + — j Pdivudx — J Ptdivudx (3.18) 


= - / p\u\ 2 dx + C / p\u\~\Vu\ 2 dx 


d f ( 1 

dt J \ 2(2 p - 


+ A)‘ 


P 2 }dx - 


1 


2(2/x + A) 


PtGdx 


Then, we will estimate the last term on the righthand side of (13.181) . 
First, it follows from (11.11) that E satisfies 


with 


( pE)t + di v(pEu) = di vF 


F = ^-V(|u| 2 ) + pu ■ Vu + Aiidivu — Pu. 


(3.19) 

(3.20) 


It follows from (13.141) and (13.191) that 

- J P t Gdx = - J (pE) t Gdx + - 2 j ( p\u\ 2 ) t Gdx 


< 


C J ( pE\u\ + |«||Vu|) \VG\dx 

1 /' 3 
-/ (div(pu)\u\ 2 G — 2 pu ■ UfG ) dx = /,. 


i= 1 


(3.21) 
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Cauchy’s and Sobolev’s inequalities yield that 

Ii+I2<v\\VG\\ 2 L 2+C(ti) f (p 2 E 2 \u\ 2 + \u\ 2 \Vu\ 2 )dx 

< ??||VG ||^2 + C{rj) J ( p\u \ 6 + p\u\ 2 + |u| 2 |Vu| 2 ) dx. 
Integration by parts and recall G = (2p + A)divu — P also gives 
h<G J p\u\ 3 \VG\dx + C J (p\u\ 2 \Vu\ + p\u\\u\) \G\dx 

<C j p\u\ 3 \VG\dx + C J (p\u\ 2 \Vu\ +p\u\\u\) (|V«| + \P\)dx 
<C(rj) J (p\u\ & + p\u\ 2 + |u| 2 |Vu| 2 ) dx + r/||VG|| 2 2 + rj J p\u\ 2 dx. 
In view of (13.161) . 


which implies 


pu = VG — pPJ x oj, 


||VG|| 2 < G||p«||| < G J p\u\ 2 dx. 


(3.22) 


(3.23) 


(3.24) 


(3.25) 


Substituting (I3.21I) - (I3.23D and (|3.25l) into (13.181) . one obtains after choosing rj suitably 
small that 

^ J &dx + J p\u\ 2 dx 

<C f (p\u\ 6 + p\u\ 2 + |u| 2 |Vu| 2 ) dx 


(3.26) 


where 


= p\S7u\ 2 + (p, + A)(divu) 2 — 2Pdivu + (2 p + A )P 2 

satisfies 

4> > Vu| 2 - C J P 2 dx. (3.27) 

Consequently, Gronwall’s inequality together with Lemma 13.21 implies (|3.17l) . This 
hnishes the proof of Lemma 13.41 

□ 


Next lemma deals with Vti. 

Lemma 3.5. Under the condition (jl.251) . it holds that for 0 < T < T*, 



j p\u\ 2 dx + J J \S7u\ 2 dxdt < C. 

(3.28) 

Proof. Make use the fact 

f = ft + div(/ (8) u) - fdivu 

(3.29) 

to obtain that 




put + pu ■ Vu + VP* + div(VP <S> u) 

= p[Aut + div(Au <8> u )] + (A + /x)[Vdi vut + div((Vdivu) <g> u)]. 


(3.30) 
















11 


(3.31) 


(3.32) 


(3.33) 


Multiplying the above equation by u and integrating over i ? 3 show that 

J p\u\ 2 dx = j Pfdivu + (u • Vu) • VPdx + p j u • [Au t + div(Au (8) u)]dx 

f 3 ^ 

+ (p + A) / ii ■ [Vdivtq + div((Vdivw) <g) u)\dx = IVj. 

J i =l 

First, recalling that 

Pt + div(Pu) + Pdivrt = 2p\D(u)\ 2 + A(divtt ) 2 
One can get after integration by parts and using the above equation that 
TV i = j P t divu + (u • Vu) • VPdx 

= - j divuAiviPu)- j PAivu-.iivvA* 

+ J (2p\D(u)\ 2 + A(divu) 2 ) diviidx + J \u • Vu) ■ V Pdx 

<j Pin- VJrfWx - / P(u • V)«dx - / 

+ C(||Vu|| £a ||V«|| La + ||Vti||i 4 ||Vu|| La ) 

< C (l + ||Vu|| 4 l4 ) + ^-||Vtt|| 2 2 

Integration by parts leads to 

JV 2 = p J Uj[dtAuj + div(uAuj)\dx 

= —p J ( diUj(diUj)t + A UjU ■ Viij) dx 

= (|v “ p " ~ * 

= —/j J (|Vu | 2 + diUjdiUjddvu — diiijdiUkdkUj — diUjdiUkdkUj) dx 

<~Y J \Vu\ 2 dx + C J |Vu| 4 Pc. 

Similarly, 

A 3 < — ^(/x + A)||divu || 2 2 + C J |Vu| 4 dx. (3.35) 

Substituting (I3.33I) - (I3.35I) into (13.311) . we obtain after choosing e suitably small that 

j t J p\u\ 2 dx + p\\Vu \\ 2 L 2 
1 j \Vu\ 4 dx + C. 


(3.34) 


(3.36) 


< C 


Note that 


||Vu|| 4 <l|Vu|| 2 ||Vu || 3 

<C||Vu|| 2 (1 + ||VG || 2 + ||Vu;|| 2 ) 

— C , ||Vu|||(l + 11 2 11 2 ) 


(3.37) 
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Substituting this estimate into (13.361) and once again recalling that 

||Vii|| 6 < ||G|| 6 + |M | 6 + C 

< C\\pu\\ L 2 +C e L\0,T), 
we conclude by Gronwall’s inequality that 

/ p\u\ 2 dx + J J \Vu\ 2 dxdt < C. 


(3.38) 


(3.39) 

□ 


Finally, the following Lemma 13.61 will deal with the higher order estimates of the 
solutions which are needed to guarantee the extension of local strong solution to be a 
global one. 

Lemma 3.6. Under the condition (jl.25|) . it holds that for 0 < T < T*, 

SU P (ll(p>-f > )ll.ff 1 rw 1 >® + l|Vw||ifi < C. (3.40) 

0 <t<T 

Proof. In view of (13.381) and (13.281) . one has 

\\Vu\\ L6 <C\\pu \\ L 2 + C <C. (3.41) 


For 2 < p < q < 6, direct calculations show that 

J t \\Vp\\ LP <C{ 1 + UVuUlcoJIIVpIIlp + C\\V 2 u\\ LP . (3.42) 

Similarly, 

^||VP||lp < c( 1 + ||Vu|| L oo)(||VP||lp + \\V 2 u\\ lp ) + C||V 2 n|| L P. (3.43) 

Applying the standard L p -estimate to (13.161) gives 

l|VG|| L 6 + ||w||i6 < C'||pfi|| L 6 < C + C '[|Vfi|| i2 , 


which shows 

\\G\\l~ < IIG'II^IIVGII^ <C + C\\Wu\\ 1 -/. 

for some (3 € (0,1). 

Applying the standard L p -estimate to (ll.lD o leads to 

||V 2 u|| L p < C (11pit11 x,p + ||VP||lp) 

<C(l + \\pu\\%\\p^-« + \\VP\\ LP ) 
<C(l + ||ph||“||Vft||^“ + ||VP|| LP ) 
<C(l + ||Vu||^ a + ||VP|| L p) 


for some a € (0,1). 
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This together with Lemma 12.41 gives 

IIVuHloc < C{ 1 + ||Vu||J/) log (e + ||V *|| L 2 + ||VP|| L ,-) + C|| Vu|| L 2 . (3.45) 

Substituting (13.441) and (13.451) into (13.42[) ~ (13.431) yields that 

/'(*) < Cg(t)f(t ) In fit), (3.46) 

where 

fit) 4 e + ||Vp|| Lff + ||VP|| l5 , g(t) 4 1 + ||Vu||£ 2 . 

It thus follows from (13.461) . (13.281) . and Gronwall’s inequality that 

sup ||(Vp,VP)|| L9 - < C, (3.47) 

0 <t<T 

which, combined with (|3.45|> and (13.281) gives directly that 

[ ||Vu||£oodi < C. (3.48) 

Jo 

Taking p = 2 in (13.421) . one can get by using (|3.48|) . (13.441) and Gronwall’s inequality 
that 

sup ||(Vp,VP)|| i2 <C, (3.49) 

0 <t<T 

which together with (|3.44l) yields that 

sup ||V 2 u || L 2 < C sup i\\pii \\ L 2 + ||Vp || L 2 + ||VP|| l2 ) < C. 

0 <t<T 0 <t<T 

This combined with (|3.47[) , (13.49[) , and (|3.17|> finishes the proof of Lemma 13.61 

□ 


Now we are in a position to prove Theorem ll.il 

Proof of Theorem 11.11 Suppose that (11.251) holds. Note that the generic con¬ 
stant C in Lemma 13.61 remains uniformly bounded for all T < T*, so the functions 

[p,u, P)(x,T*) = lim [p,u, P)(x,t) satisfy the conditions imposed on the initial data 
t->T* 

at the time t = T*. Furthermore, standard arguments yield that pii € C([0, T]; L 2 ), 
which implies 

pii{x,T*) = lim pu € L 2 . 
t^T* 


Hence, 


with 


—pAu — (p + A)Vdivw + \7P\t=r* = y/pix, T*)gix ) 


g{x) = 


p- l / 2 ix,T*)ipu)ix,T*), 

0 , 


for x € {x\p{x, T*) > 0}, 
for x € {x\p(x, T*) = 0}, 


satisfying g E L 2 due to (13.401) . Therefore, one can take ip,u, P)ix, T*) as the ini¬ 
tial data and apply Lemma 12.11 to extend the local strong solution beyond T*. This 
contradicts the assumption on T*. We thus finish the proof of Theorem 11.11 




































4 Outline of Theorem 11.2 

The main idea is quite analogous to section 3. 
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Proof. Case 7.0 = R 3 with non-vacuum far-field. 

It’s sufficient to note that 

(,p - p)t + div((/9 - p)u) + pdivu = 0. (4.1) 

Multiplying (14.11) by p — p and integrating the resulting equation over 0, we obtain 
after using Lemma 13.21 that 

(\\p-p\\hy(t)<c\\ P -p\\i 2 +c\\vu\\i 2 , 


Therefore, 


p-p€ L°°L 2 . 

(4.2) 

Similarly, one can show 


(p-p,P-P)eL oo (0,t-,L 1 nL° o ). 

(4.3) 


Then the remaining proof can be done step by step. 

Case II. n = T 3 . 

We need only to redefine the effective flux G as 

G = (2p + A)diva -P + P, (4.4) 

where 

P = i J Pdx ■ ( 4 - 5 ) 

For any q £ [l,oo) and under the condition (11.251) . G satisfies 

||G||l9 < C\\Vu\\ Lq + ||P||z,8 + C ( 4 -6) 

and 

||VG||l 9 + ||w||i9 < . (4.7) 

Then the higher regularity of the density, velocities and temperature can be obtained 
without difficulty. 

Case III. O is a bounded domain. 

Since there is no boundary condition for effective viscous flux. We will outline the 
proof of Lemma 13.41 

Motivated by [23], we decompose the velocity into two parts. It follows from Lemma 
12.41 that for any t £ [0, T], there exists a unique v(t, •) £ 77q D W satisfying 

pAv + (p, + A)Vdivu = VP, (4.8) 

which together with (12.51) yields that 


\\Vv\\ LP <C\\P\\ LP <C, for pe [2,6], ie[0,T], 


(4.9) 
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and that 


— J Ptdiwdx = — J {pEJvt ■ Vu + (p + A)dm; t divi;)dx 
= ~\j t [ (^l v?; | 2 + (/* + A)(divv) 2 ) dx. 


(4.10) 


Denoting by 

w = u — v, (4-11) 

we have w € Hq n W 2 ^, for a.e. t € [0,T], Moreover, for a.e. t € [0,T], w satisfies 

/iAui + (fi + A)Vdivu> = pu, (4-12) 

which together with the standard L 2 -estirnate for elliptic system gives 

l|Vtu|| L e + ||V 2 w;|| L 2 < C\\pu\\ L 2. (4.13) 

It follows from (13.141) and (13.191) that 

— j Ptdivwdx = — J(pE)tdivwdx + J(p\u\ 2 )tdivwdx 
< C / {pE\u\ + |it||Vit|) \V 2 w\dx 


(4.14) 


1 f ^ ^ 

- / (div(/9u)|n| 2 divr(; — 2 pu ■ utdivw ) dx = ./ ? > 


2=1 


In fact, (|4.14p has a same structure of (|3.21 D . Here divw plays a same role as G in 

mm - 

We then finish the proof of Theorem 11.21 for bounded domain by adapting a same 
procedure as Theorem 11.11 with the help of Lemma 12.41 

□ 
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